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SUMMARY 

^r^Thia  papor  doacriboa  a now  thoerotlcal  roproaontation  of  circulation  control  baaod  on  diacroto  vortex 
modolllng  tochniquoa.  Tho  application  of  diacroto  vortox  modelling  to  tho  prediction  of  circulation  control  is  initially 
preaonted  for  the  Coanda  flaw  of  a wall  Jot  around  a circular  cylindor  in  a froo  atroam.  Tho  docay  in  tho  jot 
momentum,  duo  to  vlacoua  and  entiainmont  offecta,  haa  boon  roproaontod  in  tho  modol  by  decaying  tho  atrongth  of  oach 
vortox  as  it  flows  downstream  from  tho  slot.  Tho  modol'a  application  haa  thon  boon  oxtondod  to  prodict  tho  offoct  of 
circulation  control  on  other  aha/oa  of  aorofoil.  Tho  papor  includes  tho  application  o'  tho  modol  to  an  olllptic  section  as 
an  oxamplo  of  any  aorofoil  section  that  can  bo  roproaontod  in  a potential  flow.  Tho  modol  haa  had  conaldorablo  succoss 
in  predicting  tho  lilt  produced  frf>-  t suen  bodies  duo  to  tho  jot,  but  information  about  tho  drag  roquiros  additional 
modolllng  of  th-„  surface  boundary  layer.  A similar  technique  has  boon  usod  to  roprosont  tho  boundary  layor  by  discroto 
vortices  and  doublets,  T no  modol'a  prediction  has  boon  compared  with  tho  oxporlmontal  data  of  an  unblown  circular 
cylinder. 
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Chord 

Lift  eoef.'icient 

Drag  coefficient 

Jet  momentum  coefficient 

“’art ox  distance 

Arbitrary  decay  constant 

Tangential  velocity 

Radius  of  the  cylinder 

Radial  position  of  a vortex 

Radial  position  of  an  linage  vortox 

Radial  distance  of  the  outor  edge  of  boundary  layor 

Radial  distance  of  the  outor  edge  of  displacement  thickness 

Age  of  a vortex 

Slot  thickness 

Freo  stream  velocity 

Vortex  strength 

Initial  vortex  strength 

Displacement  thickness 


1.  INTRODUCTION 

The  concept  of  circulation  control  for  high  lift  generation  was  developed  In  the  early  1960‘s  at  the  National 
Gas  Turbine  Establishment  and  haa  sinee  been  Investigated  by  a number  of  researchers.  Early  work  In  the  subject  had 
investigated  the  lift  and  drag  characteristics  of  two-dimensional  circular  and  elliptic  cross-sections  with  circulation 
control  (References  1-4).  Later  work  extended  the  application  of  circulation  control  concept  to  helicopter  rotors 
(References  5-6).  Tho  main  exigent  research  effort  la  the  X-wlng,  being  developed  for  tho  US  Navy  and  DARPA 
(References  7-9). 

Circulation  control  relies  for  lift  generation  on  the  ability  of  a wail  Jot  to  control  the  separation  point  at 
the  bluff  trailing  edge  of  an  aerofoil.  Conventional  aerofoils  feature  a sharp  trailing  edge  which  effectively  determines 
the  location  of  the  rear  stagnation  point.  The  rear  stagnation  point  of  an  aerofoil  with  a non-sharp  trailing  edge  is  not 
constrained  in  this  manner,  and  ita  location  can  be  controlled  by  a wall  jet.  Control  of  tho  aerofoil  circulation  results 
In,  what  is  known  as,  “circulation  control  by  blowing". 

The  relationship  between  the  lift  generated  by  a circulation  control  aerofoil  and  tho  wall  jet  momentum 
supplied  reflects  a complex  interaction  of  tho  aerofoil  design  parameters.  Factors  Influencing  performance  include  the 
aerofoil  shape,  leading  and  trailing  edge  design,  blowing  slot  width  and  the  number  of  slots  as  well  as  slot(s)  position. 
The  performance  of  an  arbitrary  circulation  control  aerofoil  is  determined  by  the  Interaction  of  these  design 
parameters.  If  a theoretical  model  does  not  correctly  represent  the  offocts  of  those  parameters,  It  will  fail  to  model  an 
arbitrary  circulation  control  aerofoil.  The  lack  of  a generalized  performance  prediction  method  haa  restricted 
circulation  control  development  and  applications  despite  Its  promising  performance. 

2.1  CONVENTIONAL  CIRCULATION  CONTROL  MODELLING 

Thoro  have  been  several  efforts  to  predict  tho  performance  of  an  arbitrary  circulation  control  aerofoil,  oil 
based  on  a conventional  boundary  layer  approach  to  tho  problem.  This  approach  was  developed  initially  by  Dunham 
(Reference  10)  to  predict  the  flow  about  a circular  cylinder  fitted  with  a wall  jot  exhausting  near  tangentially  to  Its 
surface.  Dunham's  method  rolled  on  the  conventional  boundary  layor  and  wall  jot  calculations  to  docido  tho  pressures  at 
tho  separation  points  on  the  uppot  and  lower  surfaces.  A solution  was  reachod  wlvon  a certain  Jet  strength  satisfied 
Thwaltes  condition  (Reference  11)  of  equal  separation  pressures  at  tho  upper  and  lower  surfaces. 
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The  model's  accuracy  of  prediction  depended  on  the  assumptions  made  regarding  the  wall  jet 
representation,  flow  entrainment  from  the  outer  flow  to  the  jet,  and  the  separation  criterion.  Dunham's  model 
satisfactorily  predicted  circulation  control  performance  trends  at  low  momentum  coefficients  (Figure  1).  However,  the 
numberical  accuracy  became  poor  at  high  momentum  coefficients  with  small  slut  widths  or  when  analysing  single  slot 
models.  This  accuracy  reflected  the  sensitivity  of  the  model  to  the  treatment  of  separation  and  entrainment.  These 
became  significant  considerations  when  the  length  of  the  surface  between  the  blowing  slot(s)  and  the  separation  point 
was  large,  either  at  high  momentum  coefficient,  or  when  only  one  slot  was  employed  at  a position  far  from  the 
separation  point. 

Similar  methods  have  been  developed  (References  12-13),  but  all  these  methods  relied  on  successful 
development  of  semi-empirical  expressions  to  treat  turbulent  boundary  layers  and  wall  jets  over  curved  surfaces  In  the 
presence  of  an  arbitrary  pressure  gradient.  Although  encouraging  results  have  been  obtained  using  this  technique,  a 
more  fundamental  approach  is  desirable.  It  is  desirable  that  such  an  approach  should  be  sufficiently  general  to  allow 
treatment  of  tho  many  applications  of  boundary  layer  control  and  Coanda  flows  in  contemporary  aerodynamics.  A new 
model  is  developed,  and  described  below,  which  is  based  on  a discrete-vortex  approximation  to  the  shear  layer  external 
to  a wall  jet  under  an  oxternal  stream. 

2.2  DISCRETE  VORTEX  MODELLING  OF  CIRCULATION  CONTROL 

Discrete-vortex  modelling  consists  of  the  representation  of  a shear  layer,  by  a number  of  discrete  vortex 
filaments  surrounded  by  an  invlscld,  irrotatlonal  flow.  Deformation  of  the  sheet  is  calculated  by  evaluating  tho  velocity 
induced  at  each  element  by  the  remaining  elements  and  satisfying  the  boundary  conditions  of  the  initial  flow.  Discrete- 
vortex  modelling  has  been  applied  to  a wide  range  of  flows  which  Include  shear  layer(s).  Rosenhead  (Reference  14) 
applied  this  technique  to  the  shear  layer  formed  between  two  uniform  streams  of  flow  with  equal  speed  In  opposite 
directions.  The  same  technique  was  used  to  represent  the  two  shear  layers  separating  a stream  from  an  exterior  stream 
flowing  with  equal  speed  In  the  opposite  direction  (Reference  13).  Gerrard  (Reference  16)  adapted  this  approach  to 
model  the  wake  shed  by  a bluff  body,  the  shear  layer  at  tho  wake  boundary  oelng  represented  by  discrete  vortices 
introduced  into  the  flow  at  successive  time  steps. 

The  application  of  this  technique  to  circulation  control  was  developed  initially  for  circular  cylinders 
(Reference  17),  then  extended  later  for  other  forms  of  aerofoils  (Reference  18).  Several  modifications  and 
developments  have  been  introduced  since,  and  the  latest  form  of  the  models  is  presented  below. 

3.1  MODEL  DESCRIPTION 

Initially,  the  flow  about  a circular  cross-section  in  a uniform  stream  Is  modelled  with  zero  circulation.  The 
circular  boundary  is  formed  by  introducing  a doublet  into  the  uniform  stream.  The  shear  layer  formed  at  the  edge  of 
the  jet  is  represented  by  a series  of  discrete  vortices.  The  shear  layer  is  formed  as  the  jet  emerges  from  the  slot,  so 
the  vortices  are  Introduced  progressively  to  simulate  the  growth  of  the  jet  os  it  leaves  the  slot  and  flows  downstream. 
The  first  vortex  is  Introduced  at  a point  above  the  circular  boundary  where  the  slot  is  located.  The  distance  of  this 
point  from  the  boundary  is  the  slot  width.  The  vortex  is  accompanied  by  an  Image  vortex  inside  the  boundary,  with  tho 
same  strength  but  opposite  sense  at  rotation.  The  position  of  the  image  vortex  is  decided  by  the  equation) 

RvRi  » R20  (1) 

where  Rv  i Radial  distance  of  the  vortex 

Rj  i Radial  distance  of  the  image  vortex 

R0  i Radius  of  the  circular  boundary 

The  reason  for  introducing  the  image  vortex  Is  to  preserve  the  shape  of  the  circular  boundary  by  ensuring 
the  absence  of  any  normal  velocity  at  the  circular  boundary.  The  external  vortex  moves  a specified  distance  step  after 
calculating  the  induced  velocities  on  this  vortex  by  all  elements  of  the  fiow  field,  including  Its  linage.  The  linage 
vortex  position  is  readjusted  according  to  equatl  n (1),  white  a new  vortex  is  Introduced  at  the  slot  lip  with  its  new 
linage.  The  process  la  repeated,  allowing  each  vortex  to  move  a distance  proportional  to  the  time  taken  by  the  initial 
vortex  to  move  over  the  fixed  distance-step,  and  also  to  the  local  velocity  of  the  vortex  in  question.  After  each  step,  a 
new  vortex/image  is  introduced  at  the  slot  Up.  The  process  leads  eventually  to  a sheet  of  discrete  vortices  as  shown  in 
Flg.2. 


3.2  JET  DECAY 

In  the  real  flow,  the  jet  momentum  decays  die  to  viscous  and  mass  entrainment  effects.  As  a segment  of 
the  jet  mass  flaw  leaves  the  slot  and  flows  downstream,  a gradual  loss  of  its  momentum  occurs  due  to  viscosity  and 
entrainment  from  the  external  flow.  This  loss  in  momentum  may  be  represented  in  the  modei  by  the  use  of  a time 
dependent  vortex  decay.  The  rate  of  loss  of  momentum  is  expected  to  be  in  proportion  to  the  initial  Jet  momentum  and 
the  external  flow  momentum.  Therefore,  an  arbitrary  form  of  decay  has  been  adopted  in  the  modei  by  reducing  the 
vortex  strength  as  it  flows  from  the  slot.  Assuming  that  the  initial  strength  of  the  vortex  is  6^),  then  its  strength  wlU 
be  reduced  to  the  value 

>(t)  - y0.exp(-  KtVr0tAl>  (2) 

where  U«i  Free  stream  velocity 
d i Distance  step 
t : Time  ago  of  the  vortex 
K i Arbitrary  constant 
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The  value (\,/d)representa  the  circulation  per  unit  length  of  the  shear  layer  when  the  jet  leaves  the  slot,  hence  It  can  be 
used  as  a measure  of  the  jet  momentum,  the  value  vm  Is  a measure  of  the  momentum  of  the  external  flow.  The  chosen 
form  of  decay  obeys  the  observed  conditions  that  the  strength  of  the  vortex  becomes  zero  as  Its  age  Increases 
indefinitely  and  that  Its  decay  Is  most  rapid,  when  the  vortex  age  is  small  and  the  shear  is  greatest. 


3.3  CALCULATION  OF  THE  LIFT  AND  JET  MOMENTUM  COEFFICIENTS  (Cl  and  Cj) 

The  lift  coefficient  (Cl)  is  calculated  by  integrating  the  normal  local  pressure  coefficients  (oJ  on  the 
circular  boundary,  determined  as  usual  from 

Si  " 1 - (u*/u-)2 

whore  u4  = local  velocity  induced  by  all  the  flow  elements. 

The  jet  mass  flow  rate  (in  ) Is  defined  as  pVjtj,  where  Vj  is  the  mean  jet  velocity  at  the  slot  exit,  tj  Is  the 
jet  thickness  at  the  slot.Vj  is  calculated  in  the  model  as  tho  surface  velocity  beneath  the  slot  lip  and  defined  as  Vs.  The 
jet  momentum  is  then  defined  as  the  product  of  rti  and  the  increment  in  surface  velocity  beneath  the  slot(AV)due  to  the 
presence  of  the  jet.  In  defining  this  jet  momentum,  a primary  consideration  is  that  the  momentum  coefficient 
predicted  should  be  zero  when  no  air  is  supplied  to  the  blowing  system.  Therefore  the  jet  momentum  coefficient  is 
calculated  as: 

Cj  - (Vs  AV  tj)/(u£  Ro) 

where  V„  is  the  total  surface  velocity  at  the  slot  and  AV  is  the  velocity  induced  at  the  surface  by  the 
vortices  only. 

4.1  MODEL  CONVERGENCE  AND  SENSITIVITY  TO  VORTEX  SPACING 

The  lift  and  momentum  coefficients  will  be  reliably  predicted  when  the  following  conditions  are  met:- 

1)  the  geometry  of  the  jet  sheet  should  stabilise  so  that  the  introduction  of  additional  vortices  have  no 

effect  on  U 

2)  the  vortex  spacing  used  to  model  the  jet  sheet  should  not  affect  tho  lift  and  momentum  coefficients. 

In  the  previous  sections  it  has  been  shown  that  the  process  of  introducing  vortices  at  the  slot  Up  eventually 
led  to  a fixed  geometry  for  the  jet  sheet  which  was  unaffected  by  tho  introduction  of  further  vortices  at  the  slot  Up, 
thus  satisfying  condition  (1)  above. 

Condition  2 was  tested  by  evaluating  Ci  and  Cj  for  different  values  of  the  vortex-distance  ratio  (d1  a d/R) 
for  appropriate  initial  vortex  strengths  determined  so  thaty  g/dLL  remained  constant.  Figures  3 and  4 show  that  the 
predicted  values  of  Cl  and  C j Is  strongly  dependent  on  the  vortex-distance  ratio  d*.  The  curves  of  both  figures  suggest 
that  there  are  unique  values  of  Cl  and  Cj  when  d*  a Q.  The  curves  do  not  continue  to  d*  a 0 because  tho  effect  that  one 
vortex  induced  on  Its  Immediate  neighbours  became  very  large  at  smell  d*  and  the  resultant  motions  of  the  sheet  are 
best  described  as  chaotic.  To  solve  this  problem,  a vortex  splitting  process  has  been  adopted,  as  explained  below. 

4.2  VORTEX  SPLITTING 

In  order  to  study  the  convergence  of  the  model's  prediction,  further  reduction  In  the  vortex  distance  is 
needed,  without  altering  the  contour  of  the  vortex  sheet.  Reduction  of  the  vortex  distance  Is  possible  by  splitting  the 
vortices  after  ties  sheet's  contour  has  been  established.  Each  vortex  in  the  sheet  Is  divided  Into  a group  of  sub-vortices. 
The  sub-vortlces  are  spread  along  de  Une  separating  this  vortex  from  the  next  one.  The  strength  of  each  sub-vortex  is 
calculated  In  proportion  with  the  general  decay  of  the  sheet  as  expressed  in  equation  (2),  provided  that  the  total 
strength  of  a group  of  sub-vortices  equals  the  strength  of  the  original  vortex  which  was  split.  This  lias  similarities  to 
Mttskew*s  work  (Ref.  19). 

The  splitting  process  allows  an  indefinite  reduction  in  the  vortex  distance  without  altering  the  contour  of 
the  vortex  sheet,  or  the  circulation. 

Figures  3 and  6 show  the  results  obtained  by  introducing  the  splitting  process  to  Ute  model.  The  figures 
show  that  both  Cl  and  Cj  continue  to  converge  towards  ultimate  values,  as  suggested  before  by  Figures  3 end  4.  A 
wider  range  of  investigation  on  different  circular  cylinders  has  concluded  that  if  the  vortex  distance  is  eventually 
reduced  to  0.1  per  cent  of  the  cylinder's  radius,  the  predicted  Cl  and  C j values  will  converge  to  within  0.3  per  cent  of 
its  ultimate  values.  The  ultimate  values  are  those  expected  if  the  vortex  distance  is  re  due  at  to  zero,  l.e.  when  the 
vortex  sheet  becomes  3 continuous  sheet  of  vorticity.  The  model,  in  this  form,  has  successfully  represented  a 
continuous  si  wet  of  vorticity  by  a discrete  vortex  model. 

34  COMPARISONS  WITH  EXPERIMENTAL  DATA 

The  mat  ha  latlcal  model  has  been  spplied  to  a wide  range  of  circular  cylinders  and  Us  predictions 
compared  with  the  cun-spending  experimental  data.  Tlve  application  has  covered  different  circular  cylinder  designs 
and  flow  conditions.  A v iirlety  of  circular  cylinder*  have  been  tested  and  reported  In  Ref.(Z).  Two  of  those  models  have 
been  chosen  to  demonstrate  the  vortex  model's  prediction  capability.  Figure  7 shows  the  comparison  for  a single  slot 
aodel  (ROMNEY  MODEL  2 OULD  4)  at  two  angle*  of  attack.  Figure  7a  represents  the  case  when  the  slot  was 
located  at  the  top  dead  centre,  with  the  wind  tunnel  speed  at  M.«  0.2.  In  Figure  7b  the  slot  was  moved  2iP  aft  from  top 
dead  centre,  while  the  tunnel  speed  was  at  s 0.34.  The  prediction  has  been  calculated  In  each  case  for  different 
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values  of  the  arbitrary  constant  (K),  defined  In  the  decay  equation  (2).  The  model  gives  an  excellent  prediction  of  the 
experimental  results  when  K = 0.16$.  Despite  the  differences  in  free  rtream  velocity  and  slot  position,  the  same  value 
of  K has  produced  equally  good  agreement  with  the  experiment  shown  in  the  second  case  (7b). 

A different  model  (ROMNEY  MODEL  2 BUILD  2)  has  been  used  for  further  comparison.  The  model 
employed  two  slots,  one  at  top  dead  centre  and  the  other  at  30°  downstream,  while  the  wind  tunnel  speeds  were  M.= 
0.2  and  0.32,  as  shown  In  Figure  (8).  The  discrete  vortex  model  has  been  applied  In  this  case  by  Introducing  vortices 
from  the  two  slots  simultaneously.  Both  sheets  have  the  same  Intensity  of  vorticlty,  the  total  jet  momentum 
coefficient  is  the  sum  of  both  jet  momentum  coefficients.  The  motion  of  each  vortex  is  decided  by  the  velocity  Induced 
on  this  vortex  by  all  vortices  in  both  sheets.  Figure  8 shows  that  the  prediction  is  consistent  with  that  produced  in  the 
single  slot  model,  since  the  value  of  K remains  virtually  unchanged.  The  comparison  of  the  discrete  vortex  model's 
prediction  with  other  circular  cylinders  has  produced  the  same  conclusion.  The  value  cf  K remains  near  constant  in  all 
cases,  regardless  of  the  model's  dimension,  slot(s)  position  on  free  stream  velocity.  A similar  conclusion  has  been 
obtained  when  the  model  is  applied  to  cross-sections  other  than  circular.  An  example  of  this  application  Is  shown 
below. 

5.2  APPLICATION  OF  THE  MODEL  TO  ELLIPTIC-SECTIONS 

The  discrete-vortex  modelling  may  be  used  to  represent  a jet  on  other  cross-sections.  The  modelling  can 
bo  applied  to  any  section  provided  that  it  can  be  represented  by  potential  flow  elements  such  as  doublets,  sinks,  sources, 
etc.  This  representation  is  essential  because  it  is  used  as  a starting  point  for  the  discrete  vortex  modelling.  However, 
an  alternative  method  to  predict  the  performance  of  a specified  cross-section  Is  to  use  the  conformal  transformation. 
If  the  required  section  can  be  produced  by  a suitable  conformal  transformation  from  a circular  section,  then  the 
velocity  distribution,  as  predicted  from  the  model  oiscussed  above,  can  be  transformed  to  predict  the  performance  of 
the  required  section.  This  method  has  been  used  to  predict  the  performance  of  an  elliptic-section.  The  Kutta- 
Zhukovsky  transformation  has  been  used  to  calculate  the  velocity  distribution  on  the  elliptic  boundary  which  allows  the 
calculation  of  Doth  lift  and  jet  momentum  coefficients  of  the  ellipse. 

A comparison  with  experimental  results  is  sliown  in  F igure  9.  The  experimental  model  is  Allcock's  elliptic 
section  10/A  (Reference  3),  the  section  has  a thickness  chord  ratio  of  7.0%,  with  the  slot  located  at  a position  of  Q.874C 
from  the  leading  edge,  C being  the  chord  of  tho  ellipse.  Figure  9 shows  a similar  result  to  that  demonstrated  for  the 
circular  section.  A value  of  K within  the  range  of  0.15-0.165  brings  the  model's  prediction  into  very  good  agreement 
with  the  experimental  data,  further  supporting  the  thought  that  K may  be  a constant  irrespective  of  section. 

6.1  BOUNDARY  LAYER  MODELLING  BY  DISCRETE- VORTICES 

The  previous  discussion  has  shown  that  the  discrete  vortox  modelling  of  the  jot's  shear  layer  successfully 
predicts  the  lift  produced  from  an  arbitrary  aerofoil  section  with  circulation  control.  The  model  in  its  present  form 
cannot  provide  a prediction  of  the  drag  produced  from  such  an  aerofoil.  An  accurate  prediction  of  the  drag  requires 
additional  modelling  of  the  surface  boundary  layer.  The  surface  boundary  layer  is  a shear  layer  which  forms  on  the 
aerofoil  boundary  duo  to  the  "non-slip"  condition  at  the  surface.  The  surface  shear  layer  la  formed  between  a tero- 
veleclty  surface  and  an  external  flow  of  non-rero  velocity,  while  the  jet  shear  layer  is  formed  between  a high  velocity 
jot  flow  and  an  external  slower  velocity  flow.  Therefore  the  surface  shear  layer  can  be  considered  aa  a sheet  of 
vorticlty  similar  to  that  of  the  jet  but  with  the  opposite  sense  of  circulation.  White  the  jet  sheet  speeds  up  the  flow  at 
the  upper  surface,  the  boundary  layer  sheet  tends  to  alow  the  flow  down.  Therefore  the  circulation  In  the  two  sheets 
will  be  opposite.  A discrete  vortex  model  of  the  boundary  layer  on  an  unblown  circular  section  lias  been  developed  prior 
to  combining  it  with  the  jet  model  in  order  to  predict  both  lift  and  drag. 

6.2  UNBLOWN  CIRCULAR  SECTION  BOUNDARY  LAYER  MODEL 

A model  of  discrete  vortices,  which  represent  the  boundary  layer,  must  satisfy  the  boundary  conditions  of 
the  flow  In  this  layer.  The  velocity  profile  Inside  the  boundary  layer  should  satisfy  the  two  basic  conditions  of  rero 
velocity  at  the  surface  and  a velocity  at  live  outer  edge  of  the  profile  equal  to  that  of  the  potential  flow.  The  presence 
of  the  vortices  should  create  a velocity  field  which,  when  combined  with  the  potential  flow  field,  leads  to  a velocity 
profile  similar  to  that  measured  In  the  boundary  layer.  Figure  10  explains  how  a pair  of  vortices  (vortex  * Image) 
create  a velocity  profile  that  otters  the  potential  flow  profile  Into  one  similar  to  that  of  a boundary  layer.  The 
magnitude  of  the  velocities  In  the  vortex  profllo  depend  on  the  strength  and  position  of  tho  vortex  relative  to  the 
circular  boundary,  as  well  as  the  angular  distance  of  the  vortex  from  the  control  point  A. 

To  calculate  the  strength  of  the  vortex,  consider  a control  area  aa  shown  in  Figure  11.  The  area  lies 
between  two  arbitrary  angular  positions  (8i.  82),  Us  Internal  edge  Is  the  circular  cylinder  surface  and  the  outer  edge 
Uee  In  the  flow  just  outside  the  boundary  layer.  The  total  circulation  bound*!  by  this  area  < Nvlcuiated  ast 

Y * | q 4s 

where  q a tie  velocity  tangential  to  the  area's  boundary 
and  ds  » element  of  distwvce  along  the  area's  boundary 

Since  q a (1  along  the  circle  surface 
*b 

1 

where  qt  1 Potential  flow  tangential  velocity 

i Normal  velocity  Inside  the  boundary  layer 
8 1 Angular  position  measured  from  the  tending  edge 
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Aaaumlng  that  the  change  in  the  normal  velocity  between  positional  and 62  is  negligible,  i.e.  4,  = a 
f«2  2 


f9  2 

It  de 

'81 


(3) 


The  circulation  calculated  in  this  way  aatlafies  the  condltiona  of  zero  aurface  velocity,  and  that  the 
velocity  outside  the  boundary  layer  is  that  of  the  potential  flow.  However,  if  equation  (3)  Is  applied  for  the  vortices  at 
positions  ‘ 3 >90®,  the  calculated  circulation  will  decrease  as  e Increases  because  of  the  reduction  in  qt . This  is  not 
physically  ceptabie,  since  the  boundary  layer  vorticlty  has  to  Increase  because  it  represents  the  momentum  destroyed 
by  the  surface.  Therefore  it  has  been  assumed  that  the  vortex  strength  remains  constant  for  the  region  8 > 90*.  The 
complete  modelling  of  the  boundary  layer  starts  by  dividing  the  layer  along  the  surface  into  (n)  control  areas  of  equal 
angular  width.  The  circulation  within  each  box  Is  calculated  according  to  equation  0)  or  according  to  the  assumption 
mentioned  above,  then  each  calculated  circulation  is  concentrated  in  a single  vortex  at  the  middle  angular  position 
between  81  and  62  (Flg.ll).  The  distance  of  each  vortex  from  the  circular  aurface  is  decided  by  applying  the  condition  of 
zero  velocity  at  (n)  control  points  on  the  circle.  The  control  points  are  chosen  at  the  intermediate  angular  positions 
between  the  vortices.  These  points  are  used  as  reference  points  for  calculating  the  velocity  profiles  across  the 
boundary  layer.  The  application  of  zero  velocity  condition  leads  to  (n)  algebraic  equations  which  can  be  solved  to 
calculate  the  exact  positions  of  the  vortices.  The  model  adopts  two  groups  of  (n)  discrete  vortices  to  represent  the 
boundary  layer  at  the  upper  and  lower  surfaces,  but  only  (h)  algebraic  equations  are  solved  because  of  the  symmetry  of 
the  flow  on  both  surfaces.  The  velocity  profile  has  been  calculated  at  each  control  point,  and  the  outer  edge  of  the 
profile  (Rfc)  is  considered  to  be  reached  when  the  calculated  velocity  equals  that  of  the  potential  flow. 


The  calculated  velocity  profile  has  been  used  to  calculate  the  displacement  thickness  at  each  control  point. 
The  displacement  thickness  is  calculated  from  the  equations: 


where  *d  * *0  * 4d 


where  dy  t distance  element  along  the  radial  direction 
: Displacement  thickness 
PF  denotes  the  original  potential  flow 
BL  denotes  the  boundary  layer  due  to  vortex  flow 


This  equation  la  driven  from  the  definition  of  the  displacement  thickneaa  as  the  distance  by  which  the 
potential  flow  streamlines  are  displaced  by  the  presence  of  the  boundary  layer. 


The  pressure  coefficient  at  this  point  is  then  defined  as 

Cp=  1.  (qtj/O)2 

where  is  the  tangential  potential  flow  velocity  at  the  (^stance  Ftp.  The  pressure  coefficient  distribution  along  the 
surface  Is  used  to  correlate  the  model's  prediction  with  the  experimental  data.  To  complete  modelling  the  flow  around 
the  circular  cylinder,  It  1s  necessary  to  model  the  separated  flow  region  behind  the  cylinder.  The  position  of  the 
separation  point  is  decided  from  the  calculated  velocity  profiles  by  using  the  definition  of  separation.  The  separation 
point  It  defined  as  the  point  where  the  surface  shear  stress  vanishes,  l.e.  when  the  tangential  velocity  gradient  at  the 
surface  becomes  zero  0.  This  gradient  hat  been  calculated  from  the  velocity  profile  at  each  control  point  until 

Its  value  reaches  zero.  The  separstion  Is  considered  to  occur  at  thia  position  and  the  separation  pressure  coefficient  is 
taken  as  the  pressure  coefficient  of  this  point.  Since  the  separated  flow  region  hat  constant  pressure,  It  la  simulated  by 
introducing  a number  of  doublets  on  the  surface  between  the  separation  points.  The  strength  of  these  doublets  ore 
decided  by  applying  the  condition  of  constant  separation  pressure  at  the  control  points  in  the  separated  area. 


The  model's  prediction  of  the  pressure  coefficient  distribution  has  been  compared  with  tcketthaeh'sdsia 
(Reference  20)  far  a supercritical  flow,  l.e.  the  turbulent  boundary  layer  la  fully  developed.  Figure  12  show*  the 
comparison  between  experiment,  potential  flow  and  model's  prediction.  The  figure  shows  how  the  number  of  vortices 
has  affected  the  predicted  separation  position  and  pressure.  It  has  been  found,  ea  suggested  by  the  figure,  that  the 
prediction  converge*  rapidly  as  tte  number  of  vortices  increase*.  If  the  number  is  increased  more  than  (0,  there  is 
Uttle  effect  on  the  predieted  position  and  pressure  at  separation.  Apart  from  the  region  between  8 • IS*  and  IQS*,  the 
prediction  has  produced  a very  good  agreement  with  the  experimental  data.  The  deficiency  In  the  prediction  in  this 
region  is  probably  dua  to  the  assumption  of  constant  vorticlty  in  the  region  of  e > 90*.  The  modal's  prediction  Is 
expected  to  Improve  if  the  vorticlty  Increases  in  this  region  rather  than  remaining  constant.  Further  investigation  it 
required  to  find  a formula  governing  the  growth  of  vorticlty  in  this  region.  The  drag  coefficient  Cq  has  been  calculated 
by  Integrating  the  predicted  preeaure  distribution,  and  compared  with  the  experimental  value.  The  predieted  pressure 
distribution  gave  a value  of  Cq  » 0.768,  white  the  corresponding  experimental  value  was  Cq  • 0.7*9,  that  la,  an  error  of 
teas  than  2%,  The  notieable  error  in  the  modal's  prediction  of  the  Cp  distribution  In  the  region  of  IQS*  > 6 > IS*  hat,  as 
expected,  a very  small  contribution  to  the  calculated  value  of  Cq. 

7.  CONCLUSIONS 


A theoretical  model  has  been  developed  to  predict  the  performance  ef  an  arbitrary  aerofoil  whose 
circulation  Is  controlled  by  blowing.  The  modal  relies  on  representing  the  existing  shear  lay  arts)  by  sheet  (s)  of  discrete 
vortices  rather  than  using  the  conventional  boundary  (shear)  layer  analysis.  Initially,  a modal  has  boon  developed  to 
predict  the  Uft  and  }et  momentum  coefficient*  of  a section  equipped  with  a slot  which  dischargee  a wall  tat 
downstream.  The  shear  layer  between  the  Jot  and  the  asternal  flow  haa  been  represented  by  a series  vertices 
emarglng  from  the  slot.  The  model  ha*  employed  a form  of  vortex  time  decay  to  simulate  the  real  Jot  decvy  due  to 
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viscous  and  entrainment  offocts.  It  has  been  assumed  that  tho  strength  of  decay  it  in  proportion  to  tho  oxtornal  flow 
momentum,  tho  jot  momentum  and  an  arbitrary  constant  of  decoy.  In  comparison  with  experimental  data,  tho  model's 
prediction  has  shown  that  tho  value  of  this  arbitrary  constant  is  approximately  fixed  dospito  tho  wido  variations  in 
model  geometry,  number  of  slots  and  oxtornal  flow  velocity.  In  order  to  predict  tho  drag  of  tho  circulation  controlled 
aerofoil,  a second  model  has  been  developed  to  represent  tho  boundary  layor  on  an  unblown  circular  aerofoil,  as  an 
initial  stop  towards  combining  the  two  modols  for  tho  prediction  of  both  lift  and  drag.  The  boundary  layor,  as  a shear 
layor,  has  bocn  represented  In  tho  modol  by  a series  of  discrete  vortices,  whilo  the  separated  flow  region  has  been 
simulated  by  a group  of  doublets.  Tho  strongth  of  vorticos/doubtets  has  boen  calculated  using  tho  boundary  conditions 
applicable  to  tho  boundary  layer  and  separation  criteria.  Tho  model's  prediction  of  tho  prossuro  distribution  around  the 
surface  has  produced  reasonable  agreement  with  tho  experimental  data  of  tho  unblown  circular  section.  The 
experimental  drag  coefficient  is  predicted  thoorotlcally,  with  an  error  of  only  2%. 

A combination  of  tho  jot  and  tho  unblown  modols  is  oxpected  to  rosult  In  a generalized  representation 
which  predicts  both  lift  and  drag  coefficients  of  an  arbitrary  circulation  controlled  aerofoil.  It  is  also  oxpocted  that  by 
combining  tho  jet  and  tho  boundary  layer  sets  of  vortices,  an  explanation  for  tho  nearly  constant  valuo  of  docay  In  tho 
jet  modol  may  bo  possible. 
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